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ABSTRACT 

We produce the most comprehensive public void catalog to date using the Sloan Digital Sky Survey 
Data Release 7 main sample out to redshift z = 0.2 and the Luminous Red Galaxy sample out to 
z = 0.44. Using a modified version of the parameter-free void finder Z0B0V, we fully take into account 
the presence of the survey boundary and masks. Our strategy for finding voids is thus appropriate 
for any survey configuration. We produce two distinct catalogs: a complete catalog including voids 
near any masks, which would be appropriate for void galaxy surveys, and a bias-free catalog of voids 
away from any masks, which is necessary for analyses that require a fair sampling of void shapes and 
alignments. Our discovered voids have effective radii from 5 to 135 h~ 1 Mpc. We discuss basic catalog 
statistics such as number counts and redshift distributions and describe some additional data products 
derived from our catalog, such as radial density profiles and projected density maps. We find that 
radial profiles of stacked voids show a qualitatively similar behavior across nearly two decades of void 
radii and throughout the full redshift range. 

Subject headings: cosmology: observations, cosmology: large-scale structure of universe, astronomical 
databases: catalogs 



1. INTRODUCTION 

The hierarchical clustering of matter in the Universe 
naturally leads to large under dense regions, called voids. 
Indeed, the presence of voids in the large-scale distribu- 
tion of galaxies was one of the ear ly predictions of cold 
dark matter cosmological theories ( Hausman et al.|1 983), 
and the discover y of voids in some of the fir st galaxy 
redshift sur veys (Gregory & Thompson 1978 Kirshner 
et al. 19811) quickly provided a rich source of interest. 



Tod ayT^alaxy surveys, such as th e Void Galaxy Sur- 



vey 



van de Weygaert et al. 20111) and the Las Com- 



panas Redshift Survey ( |Muller et al. |2000p routinely find 



and characterize both voids themselves and their con- 
tents for useful astrophysical and cos mological informa- 
tion (see |Thompson fc GregorypOllj for a review) . 

A combination of observations and simulations gives a 
coherent picture of void properties. In the cosmological 
constant plus cold dark matter (ACDM) picture of cos- 
mic evolution, voids — usually defined to have densities 
of 10-20% the cosmic mean — have characteristic radii 
of 10-40 /i _1 Mpc, with the smallest identifiabl e voids in 
the local universe hay ing radii ~ 7 h~ 1 Mpc ( jTikhonov 
|fc Karachentsev 2006). Early structure formation sim- 
ulations successfully reconstructed observ ed voids (Hoff- 
|man fc Shaham||l982[ |White et al.| | 1987[). Later studies" 
of vo ids from surveys such as IRAS ( jPlionis fc Basilakos 
|2002p , the 2dF Galaxy Redshift Survey QHoyle fc Vogeley] 
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2004), and the Sloan Digital Sky Survey ( [Pan et al. | 20lT|) 
confirmed these properties. Semi-analyt ic models ( |Ben 



son et al. 2003 Tinker fc Conroy 200 9J) and large-scale 



ab in i tio simulations ([p ubinski et al .|1993 



20051 ICeccarelli et al.||2006| |Park fc Lee 



Colberg et al. 
2007;^ 



kel 



et al. 2011) have further elucidated the evolution, inter- 



nal structure, and distribution of voids. 

Since voids are nearly empty, their dynamics are dom- 
inated by dark energy. Thus, they may pr ovide crucial 
probe s of primordi al density fluctua tions (Sahni et al. 
1994]), fifth forces (|Li fc Zhao|[2009t, and F(r) gravity 



mocfels jLi et al.||2(H2|) The Alcock-Paczynski test (JaT 
cock & Paczynski 11979ft can be applied to measurements 



of void ellipticities, directly probing the expansion his- 
tory of the universe (|Ryden|1995) [Ryden fc Melott| 1996 



Park fc Lee|[2007| |Biswas et ai"T|2010| |Lavaux fc W an- 



delt 2011). The internal structure of voids behaves as 



a universe in miniatu re, allowing for probes of the his- 
tory of dark e nergy (jGottlober et al. 2003 Goldberg 



Vogeley 2004). The ellipticity distribution of voids can 
provide insights into the growt h of structure and the cor- 
rectness of Gener al Relativity ( Sh oji fc Lee|2012| |Lavaux| 
|fc Wa ndelt 2010). Void orientation and spin statistics re- 
veal i nformation on large -scale tidal fields ( |Lee fc Park| 
2006; Platen et al. 2008). Understanding the locations 



and sizes of voids is also crucial for cosmic microwave 
background missions, since they affect th e CMB signal 
via the integrated Sachs- Wolfe effect (Tr hompson & Vish- 
niac|1987HVadas|1998HCruz et al.|20081|Gurzadyan et al. 



2009) |Granett et al.||2008 ) 
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The reliability of the above conclusions and predic- 
tions rests on the ability to robustly produce void cat- 
alogs from galaxy surveys. While void finders are well- 
studied in the context of the large - scale dark matter sim- 
ulations (e.g., |Colberg et al.j2005| |2008| ), few are applied 
directly to large-scale redshift surveys. The largest void 
catalogs previous to this work use void fin ders that rely 



on overlapping spheres of underdensities ( |Hoyle fc Voge- 



ley 2004 Pan et al. 2011). While simple to apply, this 



approach fails to capture the full geometry of the voids 
and relies on finely-tuned parameters to correctly capture 
them. Additionally, previous works ignore the presence 
of survey boundaries and masks and do not extend to 
the full redshift range of the available surveys. 

In this work, we describe our techniques for ac- 
counting for biases due to the presence of a survey 
boundary and masks. We employ these techniques 
along with a modifi ed version of the parameter- free 



void finder Z0B0V (Neyrinck| 2008 Lavaux fc Wandelt 



2011[ ) to pro duce a catalog o f voids from both the main 



sample dStrauss et al.||2002[) and luminous red galaxy 
(LRG) (Eisenstein et al. 2001) sample o f the Sloan Digi- 
tal Sk y Survey (SDSS) Data Release 7 (Abazajian et al. 
2009}. The void catalog we produce will be useful for 
many void-based astrophysical and cosmological studies, 
as already noted. This void catal og extends to h igher 
redshifts than other catalogs (e.g., |Pan et al.pOll ) and 
is the first to include not only main sample galaxies but 
also LRGs. 

We begin in Section [2] with a presentation of our selec- 
tion of data samples from the SDSS catalog. In Section [3] 
we describe our modifications to Z0B0V to handle the sur- 
vey boundary and masks as well as our process for elim- 
inating alignment biases in the void catalog. We char- 
acterize the demographics of our void catalog, including 
redshift-depedent number counts and radial density pro- 
files, in Section [4j We provide two examples of derived 
data products from the void catalog: radial profiles of 
stacked voids in Section [5] and projected density maps 
of stacked voids in Section [6j We discuss the potential 
for future applications in Section [7 and provide details 
of the layout of the public void catalog in the Appendix. 

2. DATA SAMPLES 

We identify voids in both t he SDSS main galaxy red- 
shift survey ([Strauss et al.| |2002) and the SPSS Lu- 
minous Red Galaxy (LRG) redshift survey (|Eisenstein 



et al. 



20011). We take our main galaxy sample from 
the New York Universit y Value-Added Galaxy Cata- 
log ( jBlanton et a l. 2005) which cross-matches galaxies 



from SDSS ([Abazajia n et al.||2009| ) wit h other surveys 
using impro ved photometric calibrations ( |Padmanabhan| 
et al. 2008). We draw our samples from the fulll cat- 
alog, which enforces evolution- and K-corrected magni- 
tude limits of —23 < M r < —17. This catalog extends 
to z ~ 0.3 and contains 671 , 451 g alaxies. We take the 
LRG catalog of Kazin et al. (2010), which extends from 
z = 0.16 to z = 0.47 and enforces magnitude limits of 
-23.2 < M r < -21.2. This catalog contains 105,831 
galaxies. 

The properties of a void necessarily depend on the 
galaxy distribution used to define it. We wish to have 
statistically uniform void populations, so for this work 
we choose statistically uniform, volume-limited subsets 



of the above catalog. Although our algorithm does 
not strictly require it, using volume-limited samples 
maintains similar effects of shot noise and galaxy bias 
throughout each void catalog and in principle allows 
one to make prompt comparison to predictions based on 
halo occupation models fit to galaxy correlation func- 
tions flBerlind fc Weinberg|[2002]|Zehavi et al.||2011 ). 

For the main sample, we apply evolution and K- 
corrections and compute absolute magnitudes M r assum- 
ing a co smology consistent wit h the latest WMAP 7- year 
results dKomatsu et ajL||2011| ): Q M = 0.27, ^ A = 0.73, 



and h = 0.71. Figure [T] shows the absolute r-band mag- 
nitude versus redshift for each catalog and illustrates our 
chosen samples. We choose four redshift bins. Each red- 
shift range is characterized by a typical galaxy luminos- 
ity, which we differentiate by the labels dim and bright. 
Our redshift bins are: 0.0 < z < 0.05, which we label 
diml, 0.05 < z < 0.1, labeled as dim2, 0.1 < z < 0.15, 
called bright 1, and 0.15 < z < 0.2, which we label 
bright 2. To gather as many voids as possible, we take 
each sample from redshift z = 0.0 to its given redshift 
cutoff. However, we drop from the sample any void 
whose center crosses the outer redshift limit of an in- 
terior sample. This minimizes boundary effects. For 
the LRGs, we use the tw o volume-limited s amples cal- 
culated and provided by Kazin et al.| ( |2010[ ), which we 
label here as Irgdim and Irg bright. Table [I] lists the sam- 
ple name, source catalog, limiting absolute magnitude, 
redshift bound, number of galaxies, and average galaxy 
spacing. The average galaxy spacing is (n^/V) 1 / 3 , where 
n g is the number of galaxies within each sample and V 
is the sample volume. 

3. VOIDS IN SURVEY DATA 

3.1. Coordinates 

Given a galaxy's sky latitude sky longitude </>, and 
redshift z, we transform to a hybrid coordinate system 

>cos#, 



)COS t 



where c is the speed of light and Hq is the Hubble param- 
eter at redshift z = 0. This coordinate system preserves 
relative distances; thus, we are essentially finding voids in 
redshift space. This choice is motivated by our d esire to 
apply the Alcock-Paczynski test to void shapes (Ryden 
1995) in a forthcoming work. For completeness, we will 
also make publicly available a void catalog derived from 
galaxies in real space with the same cosmology used to 
construct the volume-limited samples above. Our sam- 
ples are no longer strictly volume-limited with this choice 
of coordinates; however, the void properties that we are 
interested in should be only mildly sensitive to this in 
the narrow range of each redshift bin. Also, as we will 
discuss below the Voronoi tessellation preserves topolog- 
ical information. This means that mild distortions due 
to coordinate transformations will not destroy voids. 
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3.2. Defining voids 

We identify voids using a substantially mo dified ver- 
sion of the parameter-free void finder Z0B0V (Neyrinck 
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diml (9%) 
dim2 (23%; 
brightl (16%) 
bright2 (6%) 
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Irgdim (63%) 
Irgbright (14%) 




Red shift 



Fig. 1. — Selection of volume-limited samples. Volume-limited samples of the main sample (left) and LRGs (right). Catalog j 
are shown in black. We reduce the main sample galaxies by a factor of ten to improve clarity. We show the percentage of galaxies 
sample with the sample name. 



galaxies 
in each 



TABLE 1 

Volume-limited samples used in this work. 



Sample Name 


Catalog 


^r,max 


^-min 


^max 


Number of Galaxies 


Mean Spacing (h x Mpc) 


diml 


NYU VAGC 


-18.9 


0.0 


0.05 


63639 


3 


dim2 


NYU VAGC 


-20.4 


0.05 


0.1 


156266 


5 


brightl 


NYU VAGC 


-21.35 


0.1 


0.15 


113713 


8 


bright2 


NYU VAGC 


-22.05 


0.15 


0.2 


43340 


13 


Irgdim 


LRGs 


-21.2 


0.16 


0.36 


67567 


24 


Irgbright 


LRGs 


-21.8 


0.36 


0.44 


15212 


38 



2008 Lavaux & Wandelt 2011), which is based on a 
Voronoi tessellation of the tracer particles (in this case, 
galaxies) to reconstruct the d ensity field (jvan de Wey-| 



gaert 2007 Platen et al. 201 1[ ) followed by a watershed 



algorithm to group Voronoi cells into zones and subse- 



quently voids ( Platen et al.||2007 ). 

The algorithm works as follows. First, we build a De- 
launay tessellation of the volume from the tracer posi- 
tions in redshift coordinates. Second, we assign a density 
to each galaxy based on its Voronoi volume. Using this 
estimated density Z0B0V executes a procedure similar to 
the Watershed algorithm: the entire volume is split into 
"zones" , with each zone corresponding to the attraction 
basin of a local minimum of the density field. Z0B0V then 
assigns a core density to each zone that corresponds to 
the lowest local minimum and assembles these zones into 
voids by successively joining pairs of zones if they share 
the lowest common saddle point in the density field. This 
approach is related to the theory of per sistence brough t 
to the field of large scale structures by Sousbie| ( |2011[ ) . 
Note that there are no parameters that control the den- 
sity determination, the construction of zones, their evo- 
lution into voids, or indeed any portion of this algorithm. 

Under this framework a "void" is a collection of zones 
that share common saddle points; walls, filaments, and 
clusters naturally divide a given volume into these voids. 
The volume of the void is defined as the sum of the 
Voronoi volumes of its zones. We may record these voids 
either by their geometrical definition (e.g., a collection of 
Voronoi cells) or by the tracers contained within them. 

While all available void finders accurately identify the 
largest voids in simulations, there are many differences 



at the smallest void sizes (Colberg et al. 2008), and all 
potentially suff er from effects such as artificial bridging 
between voids ( |Park fc Lee 2009). However, methods 
based on Voronoi tessellation oner several advantages 



over competing methods (e.g., Aikio fc Mahonen 1998 
Nelson|[2009l JForero- Romero et al.||2009MGaite 



Foster 



2009r|Way et al.[|20TI| ), which typically gather voids by 
building overlapping spheres of underdensities. First, a 
Voronoi tessellation provides a natural construction of 
the local density of each particle, and hence an accu- 
rate measurement of the void volume. This construction 
is directly related to the local number of neighbors of 
a particle: more local neighbors induce smaller Voronoi 
cells and thus a higher local density, as expected. The 
use of topological criteria, like local minima and sad- 
dle points, is what mak es a void finder like Z0B0V or 
DisPerSe ( |Sousbie|2011[ ) particularly appealing for mea- 
suring geometry distortions. Since the topology is re- 
silient to metric transformations, these methods recover 
mostly the same structures even if those structures have 
been affected by redshift distortions. This resilience is 
essential for our analysis because we will be eliminating 
potential voids and arranging them in a tree based on 
density and volume characteristics. Secondly, the tessel- 
lation procedure allows voids to assume arbitrary shapes 
while still obeying an overall mean density threshold dur- 
ing the watershed phase. Finally, Z0B0V is able to work 
cleanly with any sampling of the density field — from 
sparse LRGs in surveys to dense dark matter particles in 
simulations — without fine-tuning or parameter adjust- 
ments. 

Analogous to groups and clusters, voids are naturally 
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of 



^avaux & Wandelt 



organized into a tree-like hierar chy ( Aragon-Calvo & Sza- 
lay|2012| Paranjape et al. 2012). Following the technique 



( ]201ip , we organize voids into a 
such a tree based on the natural definitions of boundaries 
and basins pro vided by Z0B0V (see Figure 1 of |Lavaux & 



| Wandelt]| 20 The tree is built from the smallest void 
to the largest. A given void accepts another void as a 
parent in the tree if all the zones of the void are present 
in the parent void and the parent void is larger than 
the considered void. This structure allows us to double 
count regions that are sampled by different void sizes at 
the same time. This is particularly useful when grouping 
void sizes into large bins. 

We impose two independent density-based selection 
criteria. First, we restrict Z0B0V to report only voids 
with mean density contrasts of —0.8 (in practice, we re- 
quire the mean density contrast of the set of Voronoi 
cells that constitute a void to be < —0.8; voids will gen- 
erally have density contrasts slightly below this thresh- 
old). This reduces the effects of sampling noise in the 
estimated size of each void ( |Schmidt et al.pOOl ). Next 
we reject any void that has overdensity greater than —0.8 
within a central region, which we take to be one quar- 
ter the effective radius (defined below). We do this to 
eliminate the effects of Poisson noise in these extremely 
underdense regions. Note that these quantities are de- 
fined relative to the mean number density of galaxies in 
the survey volume. Finally, for a given survey sample, we 
ignore any voids whose effective radius is smaller than the 
mean galaxy separation, which analy ses have shown to 



be the minimum resolvable void size flColberg et al.|2005 
Tikhonov fc Karachentsev||2006| [Platen et al.||2011[ ) 



In the following discussion, we define the center of 
each void to be the volume-weighted barycenter of all 
the galaxies contained in the void volume: 



1 



(1) 



where is the position and Voronoi volume of each 
galaxy i. In other words, this is the volume- weighted 
average position of the Voronoi cells which make up a 
void, and also the density- weighted average position of 
the galaxies. This reduces the effects of shot noise on 
the determination of the void center compared to other 
methods, such as unweighted averaging or choosing the 
most-underdense galaxy position. Also, we define the 
void radius as the effective radius, or the radius of the 
sphere which has the same volume as the Voronoi-based 
void volume: 



R 



4tt 



-V 



1/3 



(2) 



Figure [2] shows an example R = 20 h~ 1 Mpc void from 
the dim2 sample with both its galaxy members and sur- 
rounding non-members. This highlights how the Voronoi 
and watershed technique of Z0B0V stretches and deforms 
the void shape to fill the entire underdense region. The 
void is clearly buttressed by denser patches of galaxies 
identified as clusters and walls. The highly non-spherical 
shape is characteristic of voids identified with the Wa- 
tershed Transform. One consequence, which will become 
apparent below, is that a sphere of radius R is not neces- 
sarily underdense, even though the Voronoi volume itself 
has a mean density contrast of —0.8. 



3.3. Accounting for survey boundaries 

We must make a few modification to Z0B0V to account 
for survey masks, boundaries, and redshift cutoffs so that 
voids do not include volumes outside the survey extent. 
We begin with a pixelization of the survey mask using 
HEALPix QGorski et aL||2005| ). The HEALPix descrip- 
tion of the sphere provides equal-area pixels, and the 
HEALPix implementation itself provides built-in tools to 
easily determine which pixels lie on the boundary be- 
tween the survey area and any masked region. Figure [3] 
shows our pixelization of the SDSS DR7 mask and the lo- 
cation of the boundary pixels. To accurately capture the 
shape of the mask we required a resolution of n S id e = 512. 

We wish to prevent voids from extending into the 
boundary or near any masks within the survey area. To 
accomplish this we inject boundary particles along each 
boundary pixel with a density 10 _3 (/i _1 Mpc) -3 . Testing 
has shown that we require this density to ensure that ev- 
ery galaxy near the survey edge is closer to a boundary 
particle than any other galaxy. We place the bound- 
ary particles randomly within the volume defined by the 
surface area of the HEALPix pixel and the redshift ex- 
tent along the line of sight. This process and our chosen 
density results in a very thin sheath of particles that 
completely encloses the galaxy distribution. In addition, 
we place boundary particles at the minimum and max- 
imum redshift "caps" of the survey. Finally, we define 
a cubic box that completely contains the survey volume 
and distribute boundary particles evenly along the sur- 
face of that box. This last placement provides closure 
to the Voronoi tessellations of the boundary particles so 
that we do not have to directly modify this portion of 
Z0B0V. 

We assign infinite densities to the boundary particles 
— this prevents boundary particles from joining zones 
and voids and hence stops voids from expanding past 
the survey region. In addition, since the volume of any 
Voronoi cell that touches a boundary particle is by defi- 
nition arbitrary, we cannot include that cell in any con- 
structed void. Thus, after the tessellation phase we re- 
move both the boundary particles and any galaxies adja- 
cent to the boundary population (i.e., closer to a bound- 
ary particle than any other galaxy) together with their 
associated Voronoi cells. This process removed approxi- 
mately 10% of the true galaxy population. 

3.4. Handling mask-induced bias 

The above prescription produces a complete and ro- 
bust void catalog that fully accounts for survey bound- 
aries and masks: the Voronoi volumes that define each 
void lie completely within the survey volume, and the 
mean density contrast of the summed Voronoi volumes 
of each void is < 0.8. The resulting catalog is appropri- 
ate for many applications, such as studying the fractional 
volume of the Universe occupied by voids or studying the 
properties of galaxies within voids. 

However, the survey mask, boundaries, and redshift 
cutoffs introduce a subtle bias: they will preferentially 
select voids that lie along the line of sight (for bound- 
aries and masks) and perpendicular to the line of site 
(for redshift cutoffs). For example, a void parallel and 
near to the survey boundary will be relatively complete, 
but a void perpendicular to the survey boundary will ap- 
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Fig. 3. — Using HEALPix to identify boundaries. HEALPix map in 
a Mollweide projection of identified boundary zones (black) around 
and within the SDSS survey area where we inject boundary parti- 
cles. 



Fig. 2. — Example Voronoi-based void. We show the Voronoi cells that define the void in purple with galaxies in red (left). Shown is a 
void with effective radius 20 h~ 1 Mpc within a spherical region of radius 50 h~ 1 Mpc. Galaxy point sizes are proportional to their distance 
from the point of view. Galaxies interior to the void are shaded dark red. On the right is a 5 h~ 1 Mpc thick slice through the same void, 
showing exterior galaxies in black and interior galaxies in red. The orientation of the void is different between the panels to highlight 
different aspects of the structure. 

We begin with a discussion of the identified void lo- 
cations and distributions within the four volume- limited 
samples. Figures [4j |5j and [6] shows the barycenters of 
all voids in each sample overlaid on galaxy positions. To 
clarify plotting, we have rotated all galaxy and void po- 
sitions about the x-axis so that they lie on the x-y plane 
(i.e., a rotational projection). We only show galaxies and 
voids within a 25 degree opening angle. We plot the void 
positions with different colors to distinguish their sizes. 

We immediately notice that voids naturally avoid the 
edges of the survey volume. Voids that intersect the 
edges tend to be smaller, and we remove from the catalog 
any voids smaller than the mean galaxy separation. This 
is especially evident in the Irgbright subsample: since the 
mean galaxy separation is so large compared to the red- 
shift extent of the subsample, we only find the few voids 
near the median redshift. Despite this natural edge- 
avoidance, as we will see below there are still voids that 
need to be removed to produce a fully bias-free sample. 
We do not find any voids in any sample in the southern 
sky; the survey slices are so thin that they truncate any 
identified void there. 

The surviving voids distribute as expected: we clearly 
see smaller voids "hugging" the edges of filaments and 
sheets, while larger voids inhabit the more expansive 
empty regions of the survey. As we move to higher- 
redshift and sparser samples, we see ever-larger void 
sizes, although for each sample the same general distribu- 
tion of smaller and larger voids applies. In all samples, 
we tend to find more small voids near the boundaries, 
where larger voids have been truncated. 

Figure [7] shows the size distribution of voids in each 
of our samples. While we see the same number of large 
voi ds in the diml and dim2 samples as in the study 
of Pan et al. (2011), which used a single volume-limited 



pear as a smaller, truncated void, which will tend to be 
below the resolution threshold and hence expunged from 
the catalog. Thus, analyses that rely on the shapes of 
voids, such as the ellipticity probability distribution or 
the Alcock-Paczynski test, will be strongly affected be- 
cause they will see more voids aligned with the mask and 
boundary than perpendicular to it. 

To eliminate this bias, we produce a culled sample of 
voids that we refer to as the "central" sample in later 
discussions. To produce this sample we remove any void 
which, when rotated in any direction about its barycen- 
ter, intersects a boundary galaxy. Operationally we per- 
form this by taking the distance from the void barycen- 
ter to the furthermost particle and comparing that to 
the distance to the nearest boundary particle. This en- 
sures that we have a fair distribution of void shapes and 
alignments within the survey volume. 

4. VOID DEMOGRAPHICS 
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Sample: dim1 

160 i 1 1 1 




1 1 1 1 1 1 1 

-150 -100 -50 50 100 150 

x [h~ 1 Mpc] 




R - 5-1 5 h'] Mpc • R - 25-45 h~ 1 Mpc • 
R - 15-25 h" 1 Mpc • 

Fig. 4. — Void distributions throughout each dim sample. We show the spatial distribution of all voids for the diral and dim2 volume- 
limited samples. We plot galaxies in grey and void barycenters in colors depending on the size as indicated in the plots. We compact the 
galaxy and void locations by rotating their positions about the x-axis to lie within the x-y plane, and only show galaxies and voids within 
a 25 degree opening angle through the survey volume. 
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Sample: brightl 
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Fig. 5. — Void distributions throughout each bright sample. We show the spatial distribution of all voids for the brightl and bright2 
volume-limited samples. See Figure p] for a plot description. 
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Sample: Irgdim 
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Fig. 6. — Void distributions throughout each lrg sample. We show the spatial distribution of all voids for the Irgdim and Irgbright 
volume-limited samples. See Figure [4] for a plot description. 
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Fig. 7. — Distribution of void sizes. We plot histograms of void 
radii, colored by sample. Solid lines are from the all- void sample, 
while dashed lines are central voids. 

sample out to z = 0.1, we see many smaller voids because 
we are using an extra low-redshift subsample. Each suc- 
cessively higher-redshift sample produces larger voids. 
Our smallest voids, in the diml sample, have effective 
radii ~ 5 /i _1 Mpc, while the very largest in the Irgdim 
reach R = 135 /i _1 Mpc. This is caused by several fac- 
tors. First, brighter galaxies are more strongly biased, 
which can lead to larger voids at higher redshifts. Sec- 
ond, shot noise is larger in a sparser sample, and some 
void sizes may be enhanced by including regions that are 
underdense because of sampling effects, but 



Bos et al.| 

( |2012j ) find that this is a small effect. Third, at lower 
redshifts the ratio of survey surface area to volume is 
much higher; this truncates any potential large void at 
low redshift. Finally, there is a simple volume effect: if 
a void of a given size has a given probability to occupy a 
region of space due to cosmic variance, we require a cer- 
tain minimum sample volume to discover it. While there 
is some overlap of void sizes between adjacent samples, 
the multiple samples allow us to extract a large range of 
void sizes. 

By removing voids that can intersect the mask after 
rotation, we mostly affect the smaller voids of each sam- 
ple. The mask already truncates larger voids so that they 
tend to appear as smaller voids. Thus, the remaining 
large voids are unaffected during the production of the 
central catalog. This is especially true for the Irgbright 
sample, where the catalog contains only the very largest 
voids, which all survive in the central catalog. The de- 
biasing procedure removes 30 — 50% of the smaller voids 
in each sample. 

We show binned void number counts as a function of 
redshift in Figure [8j We co unt roughly half as many 
voids below redshift 0.1 as did Pan et al. (2011). This is 



not surprising since we have a strict minimum void size 
cutoff criterion that eliminates the many small spurious 
voids near the mask edges. 

Within each volume-limited sample with redshift z < 
0.2, the number of voids grows strongly with redshift. 
Since the volume-limited samples have a fixed number 
density of galaxies throughout the redshift range, we ex- 
pect to be able to count all the smaller — and more com- 
mon — voids as we probe to higher redshifts within the 
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Fig. 8. — Redshift- dependent number counts. We show binned 
number counts as a function of redshift, colored by sample. Num- 
ber counts of all void sizes as a function of redshift for each sample 
show that small voids are much more common in the data. Solid 
lines are from the all-void catalog; dashed lines are central voids. 

samples. This trend is not as strong for the LRG-based 
samples, since we can only identify very large voids. For 
the Irgbright sample, our voids are confined to a very 
narrow redshift range centered on z ~ 0.4. Within each 
sample the number density per unit volume remains rel- 
atively constant. However, between slices the number 
density drops as we lose the smallest voids due to the 
reduced galaxy number densities. 

We see the strongest effect of cutting out biased voids 
at the redshift caps of each sample: we remove 30 — 
40% of the voids at the highest redshift ranges of each 
sample and only a small fraction of the voids within the 
bulk of the sample volume. The surface area of each 
redshift boundary is much larger than the line-of-sight 
boundaries; thus, we expect to remove more truncated 
voids here. 

5. RADIAL PROFILES 

We now move to another example data product of our 
catalog: radial profiles. We build radial profiles of the 
mean density in thin spherical shells of several stacks 
of voids. Each stack contains voids with radii within a 
5 h~ 1 Mpc range. To construct the stacks, we reposition 
each void such that their bary centers coincide and rotate 
them so that their directions along the line of sight share 
a common vector. We show these profiles in Figure [9j In 
this figure, we show stacks from the all-void catalog in 
solid lines and stacks from central voids with dotted lines. 
We also show the number of voids in the given stack. To 
construct these, we do not rescale the voids, since the 
radial bin widths are small compared to the sizes of the 
voids, and rescaling voids requires difficult calibration 
of the normalization. We have normalized each density 
profile to the mean number densities of galaxies in the 
sample, n. 

These stacked void profiles clearly show a qualitatively 
similar behavior across all void sizes: an extremely un- 
derdense center (by construction), a steep wall, a large 
compensation at the wall, and a gradual decline to the 
mean density. These radial profiles reflect the qualita- 
tively unive rsal behavior of densities i nside the shell re- 
gion seen in Lavaux & Wandelt (2011). The steep walls 
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Fig. 9. — Qualitatively similar radial profiles of stacked voids. Radial profiles (i.e., mean density in thin spherical shells) of stacked 
voids for various void sizes versus the mean void radius in each stack indicate steep profiles and large compensations across all void sizes. 
Note that since different samples probe different redshift ranges, these plots also give the redshift evolution of similarly-sized voids (albeit 
identified in different types of galaxies). We indicate the number of voids in each stack with the sample name in each figure. Solid lines 
are from the all-void catalog; dashed lines are central voids. 



observed here are expected in line with other observa- 
tions: we are building profiles based on a relatively sparse 
sampling of the underdensity, so we expect voids from ob- 
servations to have sharply- defined edges and large gra- 
dients at the walls ([Fillmore fc Goldreich 1984] |Benson| 
|et al.||2QQ3| jFurlanetto & Piran||2UU6j ). Indeed, our pro- 
files are similar to those found in previous works ba sed on 
earlier observations and mock galaxy catalogs (e.g. Hoyle 



earlier observations and mock galaxy catalogs (e.g. Moyle 
| fc Vogeley l [2QQ4l |Padilla et al. | 2005j) . The density profile 



shown in |Pan et al.| ( |201l| ) does not show the characteris- 
tic overdensity because their voids are based on overlap- 
ping spherical underdensities and they average together 
and rescale voids of all sizes. Simulations, which probe 
the underlying dark matter distribution with high sam- 
pling density, produce voids with shallower walls (e.g., 
|Colberg et aJL|2QQ5j [Lavaux fc Wandert||2011| ). 

In some cases the density profile drops below the mean 
density at large distances, especially for the largest voids 
in the stack. This is especially apparent in the Irgbright 
sample, which never reaches the mean density (see the 
lower panel of Figure |6| . While we have constructed our 



catalog such that the interiors of voids do not include 
any regions outside the survey, at large radial distances 
from the void center we eventually encounter the mask, 
where the spherically- averaged mean density drops pre- 
cipitously. However, this is only a minor issue affecting 
the largest radii for all the other stacks. 

Even though we construct our voids to have mean over- 
densities of —0.8, the radial profiles reach the mean den- 
sity at much smaller radii. Since voids are ellipt ical, with 
a mean ratio of major to minor axis of two (Lavaux & 
|Wandelt]|201Q[ ), we expect spherical profiles to reach the 
void wall along the minor axis first, producing a steeper 
profile. Also, since we are stacking voids of different sizes 
without rescaling, the smaller voids will add to the den- 
sity at smaller radii, again producing steeper profiles. 

The effects of shot noise are apparent in stacks 
that contain few voids. For example, the 50 — 
55 h~ 1 Mpc stack of bright2 voids contains only five 
voids, and a single void with an excess of galaxies near 
the center drives up the stacked density profile in that 
region. 
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6. PROJECTIONS 

Our fourth and final example derived data product is 
projections of stacked voids, which are useful — after 
appropriate rescaling onto the sky — for analyses such 
as probing the ISW effect in CMB observations (e.g., 
Thompson fc Vishniac||1987| |Granett et al.| |2008) and 
reconstruction of the undistorted void shape. We show 
in Figure [lO] the projected density maps of stacked voids 
for various bin sizes. To construct these we assume a 
fiat sky approximation within the area of the projected 
stacked void. In cases where multiple samples produce 
voids within the same radial bin, we have combined their 
voids into a single projection. We take all galaxies within 
a spherical region with radius twice that of the void ef- 
fective radius. 

As with the radial profiles, we see a qualitatively sim- 
ilar distribution that roughly scales with void size: a 
minimum-density core, a strongly-defined inner wall, and 
a gradual decrease in density. The signal quickly de- 
grades as we move to lower-density samples: the voids 
structure is barely visible in the 50 — 55 /i _1 Mpc stack 
and altogether lost to Poisson noise in the 90 — 
95 h~ 1 Mpc stack. However, these projections reveal 
structures not apparent in the azimuthally-averages ra- 
dial profiles. Distortions or elongations in the wall struc- 
ture will simply appear as wider compensation regions 
in the profiles but will be immediately noticeable here. 
Also, these projections effectively remove any cosmology 
dependence of the void shape, since we have projected 
them along the redshift direction. 

7. CONCLUSIONS 

We have modified the parameter-free void finding al- 
gorithm Z0B0V to account for the survey boundaries and 
internal masks in observational data sets. This prevents 
voids from growing past the survey boundary or into any 
internal masks. Thus our approach is more generally ap- 
plicable to any given survey and mask. To demonstrate 
our technique we have constructed the first public void 
catalog using the full extent of the SDSS DR7 spectro- 
scopic survey, including the LRGs. We combined multi- 
ple volume-limited samples of the SDSS galaxy catalogs 
to maximize the number of discovered voids. We have 
produced two catalogs: one catalog that includes all dis- 



covered voids, including truncated voids near the survey 
boundaries, and a central catalog, which removes voids 
with questionable shapes and alignments. 

The general statistics of our void catalog, such as num- 
ber counts as a function of redshift and size distributions, 
broadly agree with — but significantly ex tend — both 
past analyses of observational data (e.g., Muller et al. 



[20001 |van de Weygaert et al.|2011|rPan et al.|20111|Patiri 

s from simulations (e.g., Dubinski 



|et al.||2012D and results 



|et al.|l993[|Park fc Lee|2007| ). In addition, radial profiles 
and projections of stacked voids show a qualitatively sim- 
ilar shape across the entire sample and agree well with 
previous efforts. 

Our catalogs are useful for many pursuits, including 
studies of the ellipticity distribution of voids, correla- 
tions of void positions with CMB fluctuations, Alcock- 
Paczynski tests using the shapes of voids in redshift 
space, and studies of the properties of galaxies within 
voids. We have constructed useful data sets to enable 
these studies, such as catalogs of void galaxies, void 
stacks of various radial sizes, and two-dimensional pro- 
jections of void densities. We have constructed these 
data sets using both all discovered voids and a central 
void catalog free from survey edge effects. 

We are making our catalogs and data products publicly 
available at http://www.cosmicvoids.net. 
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APPENDIX 
LAYOUT OF THE CATALOG 



We provide the catalog as a single downloadable gzip-archived file at http://www. cosmi cvoids . net] Most files are 
in human-readable ASCII format, with the exception of the raw Z0B0V- generated catalog in binary and the projections 
which are NumPy array files. 

There are two top-level directories within the catalog: figures and sdss_dr7. sdss_dr7 contains the void catalog 
and subsequent analysis as presented in this paper. The figures directory contains data files necessary to reproduce 
the figures found in this paper. 

We divide the catalog into directories based on each of our six volume-limited subsamples. Each directory contains 
the complete ZOBOV void catalog (i.e., without any size or density cuts). For each sample there are four files which 
describe the voids and their member zones and galaxies: one text file which lists the voids, two binary files which 
link voids to zones and zones to particle members, and finally a binary file which contains the galaxy positions of 
the volume-limited sample. Note that there are two void description files, corresponding to the all and central void 
catalogs. To aid in parsing these files, we provide a small exampling catalog reading routine in C, dumpVoidParticles. 

Each sample directory also contains two subdirectories: all and central, which only include voids after size and 
density cuts. Each of these subdirectories contains two text files: centers.txt and sky_positions.txt, which 
describe the coordinates and properties of each void. 

The figures directory contains data files which can be used to reproduce various figures in this paper. The file 
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Fig. 10. — Projected distributions of stacked voids. Here we show two-dimensional projections of stacked voids in various radial size bins. 
In cases where multiple samples contribute to a stack, we have combined the data. 

names indicate the figure and whether they are from the central or all catalog. 

An alternate catalog, called lcdm, is available and contains only a raw Z0B0V void catalog where galaxy positions 
are converted to real space using a ACDM cosmology. 



REFERENCES 



Abazajian, K. N. et al. 2009, ApJS, 182, 543 
Aikio, J. & Mahonen, P. 1998, ApJ, 497, 534 
Alcock, C. & Paczynski, B. 1979, Nature, 281, 358 
Aragon-Calvo, M. A. & Szalay, A. S. 2012, ArXiv e-prints: 
1203.0248 

Benson, A. J., Hoyle, F., Torres, F., & Vogeley, M. S. 2003, 

MNRAS, 340, 160 
Berlind, A. A. & Weinberg, D. H. 2002, ApJ, 575, 587 
Biswas, R., Alizadeh, E., & Wandelt, B. 2010, Phys. Rev. D, 82 
Blanton, M. R. et al. 2005, The Astronomical Journal, 129, 2562 
Bos, E. G. P., van de Weygaert, R., Dolag, K., & Pettorino, V. 

2012, ArXiv e-prints: 1205.4238 
Ceccarelli, L., Padilla, N. D., Valotto, C, & Lambas, D. G. 2006, 

MNRAS, 373, 1440 
Colberg, J. M., Sheth, R. K., Diaferio, A., Gao, L., & Yoshida, N. 

2005, MNRAS, 360, 216 
Colberg, J. M. et al. 2008, MNRAS, 387, 933 

Cruz, M., Martmez-Gonzalez, E., Vielva, P., Diego, J. M., Hobson, 
M., & Turok, N. 2008, MNRAS, 390, 913 



Dubinski, J., da Costa, L. N., Goldwirth, D. S., Lecar, M., &; Piran, 

T. 1993, ApJ, 410, 458 
Eisenstein, D. J. et al. 2001, AJ, 122, 2267 
Fillmore, J. A. & Goldreich, P. 1984, ApJ, 281, 9 
Forero- Romero, J. E., Hoffman, Y., Gottlober, S., Klypin, A., & 

Yepes, G. 2009, MNRAS, 396, 1815 
Foster, C. & Nelson, L. A. 2009, ApJ, 699, 1252 
Furlanetto, S. R. & Piran, T. 2006, MNRAS, 366, 467 
Gaite, J. 2009, JCAP, 2009, 004 

Goldberg, D. M. & Vogeley, M. S. 2004, ApJ, 605, 1 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, 

F. K., Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759 
Gottlober, S., Lokas, E. L., Klypin, A., & Hoffman, Y. 2003, 

MNRAS, 344, 715 
Granett, B. R., Neyrinck, M. C, & Szapudi, I. 2008, ApJ, 683, L99 
Gregory, S. A. & Thompson, L. A. 1978, ApJ, 222, 784 
Gurzadyan, V. G., Kashin, A. L., Khachatryan, H. G., Kocharyan, 

A. A., Poghosian, E., Vetrugno, D., & Yegorian, G. 2009, A&A, 

506, L37 

Hausman, M. A., Olson, D. W., & Roth, B. D. 1983, ApJ, 270, 351 



13 



Hoffman, Y. & Shaham, J. 1982, ApJ, 262, L23 
Hoyle, F. & Vogeley, M. S. 2004, ApJ, 607, 751 
Kazin, E. A. et al. 2010, ApJ, 710, 1444 

Kirshner, R. P., Oemler, A., Schechter, P. L., Sz Shectman, S. A. 

1981, ApJ, 248, L57 
Komatsu, E. et al. 2011, ApJS, 192, 18 

Kreckel, K., Ryan Joung, M., & Cen, R. 2011, ApJ, 735, 132 
Lavaux, G. Sz Wandelt, B. D. 2010, MNRAS, 403, 1392 
— . 2011, eprint arXiv:1110.0345 
Lee, J. Sz Park, D. 2006, ApJ, 652, 1 

Li, B., Zhao, G.-B., & Koyama, K. 2012, MNRAS, 421, 3481 
Li, B. & Zhao, H. 2009, Phys. Rev. D, 80 

Muller, V., Arbabi-Bidgoli, S., Einasto, J., & Tucker, D. 2000, 

MNRAS, 318, 280 
Neyrinck, M. C. 2008, MNRAS, 386, 2101 

Padilla, N. D., Ceccarelli, L., & Lambas, D. G. 2005, MNRAS, 363, 
977 

Padmanabhan, N. et al. 2008, ApJ, 674, 1217 

Pan, C., Vogeley, M. S., Hoyle, F., Choi, Y.-y., Sz Park, C. 2011, 
000 

Paranjape, A., Lam, T. Y., & Sheth, R. K. 2012, MNRAS, 420, 
1648 

Park, D. Sz Lee, J. 2007, ApJ, 665, 96 
— . 2009, MNRAS, 400, 1105 

Patiri, S. G., Betancort-Rijo, J., Sz Prada, F. 2012, ArXiv e-prints 
1203.6407 

Platen, E., van de Weygaert, R., Sz Jones, B. J. T. 2007, MNRAS, 
380, 551 



— . 2008, MNRAS, 387, 128 

Platen, E., van de Weygaert, R., Jones, B. J. T., Vegter, G., Sz 

Calvo, M. A. A. 2011, MNRAS, 416, 2494 
Plionis, M. Sz Basilakos, S. 2002, MNRAS, 330, 399 
Ryden, B. S. 1995, ApJ, 452, 25 
Ryden, B. S. Sz Melott, A. L. 1996, ApJ, 470, 160 
Sahni, V., Sathyaprakah, B. S., & Shandarin, S. F. 1994, ApJ, 431, 

20 

Schmidt, J. D., Ryden, B. S., & Melott, A. L. 2001, ApJ, 546, 609 

Shoji, M. & Lee, J. 2012, ArXiv e-prints 

Sousbie, T. 2011, MNRAS, 414, 350 

Strauss, M. A. et al. 2002, AJ, 124, 1810 

Thompson, K. L. & Vishniac, E. T. 1987, ApJ, 313, 517 

Thompson, L. A. Sz Gregory, S. A. 2011 

Tikhonov, A. V. & Karachentsev, I. D. 2006, ApJ, 653, 969 

Tinker, J. L. Sz Conroy, C. 2009, ApJ, 691, 633 

Vadas, S. L. 1998, MNRAS, 299, 285 

van de Weygaert, R. 2007, ArXiv e-prints: 0707.2877 

van de Weygaert, R. et al. 2011, Environment and the Formation 

of Galaxies: 30 years later, ed. I. Ferreras Sz A. Pasquali, 

Astrophysics and Space Science Proceedings (Berlin, Heidelberg: 

Springer Berlin Heidelberg) 
Way, M. J., Gazis, P. R., Sz Scargle, J. D. 2011, ApJ, 727, 48 
White, S. D. M., Frenk, C. S., Davis, M., Sz Efstathiou, G. 1987, 

ApJ, 313, 505 
Zehavi, I. et al. 2011, ApJ, 736, 59 



